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Abstract 

A generic search for anomalous production of events with at least three charged leptons is pre- 
sented. The search uses a pp-collision data sample at a center-of-mass energy of yfs = 7 TeV cor- 
responding to 4.6 fb _1 of integrated luminosity collected in 201 1 by the ATLAS detector at the CERN 
Large Hadron Collider. Events are required to contain at least two electrons or muons, while the third 
lepton may either be an additional electron or muon, or a hadronically decaying tau lepton. Events are 
categorized by the presence or absence of a reconstructed tau-lepton or Z-boson candidate decay- 
ing to leptons. No significant excess above backgrounds expected from Standard Model processes is 
observed. Results are presented as upper limits on event yields from non-Standard-Model processes 
producing at least three prompt, isolated leptons, given as functions of lower bounds on several kine- 
matic variables. Fiducial efficiencies for model testing are also provided. The use of the results 
is illustrated by setting upper limits on the production of doubly-charged Higgs bosons decaying to 
same-sign lepton pairs. 
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I. INTRODUCTION 

Events with more than two energetic, prompt, and iso- 
lated charged leptons are rarely produced at hadron col- 
liders. Such events offer a clean probe of electroweak 
processes at high center-of-mass energies, and their pro- 
duction at enhanced rates above Standard Model pre- 
dictions would constitute evidence for new phenomena. 
Models predicting events with multiple leptons in the fi- 
nal state include excited neutrino models 0, , fourth- 
generation quark models @ , the Zee-Babu neutrino mass 
model supersymmetry (7; 15 k and models with 

doubly-charged Hi ggs bosons [16j, including Higgs 
triplet models [Til Il9j|. The production of such multi- 
lcpton events in the Standard Model is dominated by 
WZ and ZZ production, where both bosons decay lep- 
tonically. Smaller contributions come from events with 
top-quark pairs produced in association with a W or Z 
boson, and from triboson production. Isolated but non- 
prompt lepton candidates misidentified as prompt arise 
in Drell-Yan events produced in association with a pho- 
ton that converts in the detector and is reconstructed as 
an electron. Prompt but non-isolated leptons misiden- 
tified as isolated can arise from Dalitz decays [20|, |21| . 
Additional non-prompt, non-isolated leptons arise from 
heavy-flavor decays and from mesons that decay in flight. 
Fake leptons can arise from hadrons that satisfy the lep- 
ton identification criteria. 

This paper presents a search for the anomalous pro- 
duction of events with at least three charged leptons in 
the final state. The search uses a data set collected in 
2011 by the ATLAS detector at the CERN Large Hadron 
Collider corresponding to 4.6 fb _1 of pp collisions at a 
center-of-mass energy of i/s = 7 TeV. Events are re- 
quired to have at least two isolated electrons or muons, 
or one of each, while the third lepton may be either an 
additional electron or muon or a hadronically decaying 
tau lepton (r ha d)- 



Selected events are grouped into four categories by the 
presence or absence of a Th a( j candidate and by the pres- 
ence or absence of a combination of leptons consistent 
with a Z-boson decay. The search is carried out sep- 
arately in each category by inspecting several variables 
of interest. The results of the search are presented as 
model-independent limits. Efficiencies for selecting lep- 
tons within the fiducial volume are also presented in order 
to aid the interpretation of the results in the context of 
specific models of new phenomena. 



Related searches for new phenomena in events with 
multilepton final states have not shown any significant 
deviation from Standard Model expectations. The CMS 
Collaboration has conducted a search similar to the one 
presented here using 4.98 fb _1 of 7 TeV data [Hj|. The 
ATLAS Collaboration has performed a search for super- 
symmetry in final states with three leptons (23|, as have 
experiments at the Tevatron [24|, HH ■ The search pre- 
sented here complements the previous searches by pro- 
viding limits outside of the context of a specific model of 
new phenomena. 



This paper is organized as follows: the ATLAS detec- 
tor is described in Section UH followed by a description of 
the samples and event selection in Sections Mil and IIV1 
respectively. The categorization of events and definition 
of signal regions is presented in Section [V] The back- 
ground estimation techniques and the results of the appli- 
cation of those techniques in control regions are described 
in Section IVII Systematic uncertainties are discussed in 
Section IVII1 The results of the search are presented in 
Section I Villi Fiducial efficiencies for model testing are 
provided in Section [IXI and are used to set upper limits 
on the pair-production of doubly-charged Higgs bosons. 
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II. THE ATLAS DETECTOR 

The ATLAS experiment [26[ is a multipurpose parti- 
cle physics detector with a forward-backward symmet- 
ric cylindrical geometry and nearly 4ir coverage in solid 
angle 27]. The inner tracking detector covers the pseu- 
dorapidity range |?y|< 2.5, and consists of a silicon pixel 
detector, a silicon microstrip detector (SCT), and, for 
\r)\< 2.0, a straw tube transition radiation tracker. The 
inner detector is surrounded by a thin superconducting 
solenoid providing a 2 T magnetic field. The calorime- 
ter system covers the pseudorapidity range \r/\ < 4.9. 
Within the region |?7| < 3.2, electromagnetic calorimetry 
is provided by barrel and end-cap high-granularity lead 
liquid-argon (LAr) electromagnetic calorimeters, with 
an additional thin LAr presampler covering < 1.8, 
to correct for energy loss in material upstream of the 
calorimeters. Hadronic calorimetry is provided by the 
steel/scintillating-tile calorimeter, segmented into three 
barrel structures within \rj\ < 1.7, and two copper/LAr 
hadronic endcap calorimeters. The solid angle coverage is 
completed with forward copper/LAr and tungsten/LAr 
calorimeter modules optimised for electromagnetic and 
hadronic measurements respectively. The muon spec- 
trometer surrounds the calorimeters. It consists of three 
large air-core superconducting toroid systems with eight 
coils each and stations of precision tracking and trigger 
chambers providing accurate muon tracking for \rj\ < 2.7. 
A three-level trigger system 28] is used to select events 
for further analysis offline. 



III. MONTE CARLO SIMULATION AND DATA 
SETS 

Monte Carlo (MC) simulation samples are used to 
estimate backgrounds from events with three prompt 
leptons. The ATLAS detector is simulated using 
GEANT4 Q , and simulated events are reconstructed us- 
ing the same software as that used for collision data. 
Small post-reconstruction corrections are applied to ac- 
count for differences in efficiency, momentum resolution 
and scale, and energy resolution and scale between data 
and simulation [3(1 l3ll] . 

The largest Standard Model backgrounds with at least 
three prompt leptons are W Z and ZZ production where 
the bosons decay leptonically. These processes are mod- 
eled with Sherpa 1.4.1 32]. These samples include the 
case where the Z boson (or 7*) is off-shell, and the 
7* has an invariant mass above twice the muon (tau) 
mass for 7* — > /i/i (7* — » tt), and above 100 MeV 
for 7* — > ee. Diagrams where a 7* is produced as 
radiation from a final-state lepton and decays to addi- 
tional leptons, i.e. W — > i*v — > £j*v — > li'i'v and 
Z -> it -> Ubf iii'i', where I and i! need not have 
the same flavor, are also included. The leading-order 
predictions from Sherpa are cross-checked with next- 
to-leading-order calculations from POWHEG-box 1.0 [11]. 



Diagrams including a Standard Model Higgs boson have 
negligible contributions in all signal regions under study. 

The production of tt + W/Z processes (also denoted 
tt + V) is simulated with MadGraph 5.1.3.28 [H for 
the matrix element and Pythia 6.425 35] for the parton 
shower and fragmentation. Corrections to the normaliza- 
tion from higher-order effects for these samples are 20% 
for tt + W {36] and 30% for tt + Z 37]. Leptons from 
Drell-Yan processes produced in association with a pho- 
ton that converts in the detector (denoted Z + 7 in the 
following) are modeled with Pythia. Additional sam- 
ples are used to model dilepton backgrounds for control 
regions with fewer than three leptons. Events from tt 
production are simulated with MC@NLO 4.01 [38l |. with 
herwig 6.520 [39| for theparton shower and fragmen- 
tation, and jimmy 4.31 [40| for the underlying event. 
Events from W+jets and W +7 production are simulated 
with ALPGEN 2.13 [41] for the matrix element, herwig 
for the parton shower and fragmentation, and JIMMY for 
the underlying event. 

Simulated samples of pair-produced doubly-charged 
Higgs bosons [H, [lj], Ell are used to illustrate the re- 
sults of this search in the context of a specific scenario. 
The doubly-charged Higgs bosons decay to pairs of same- 
sign leptons, producing up to four energetic, prompt, 
isolated charged leptons in the final state. The doubly- 
charged Higgs bosons are simulated with masses ranging 
from 100 GeV to 500 GeV. A sample of pair-produced 
fourth-generation down-type quarks [3] is also consid- 
ered when estimating fiducial efficiencies and potential 
contributions from non-Standard-Model processes. In 
this model, the heavy quarks decay to top quarks and 
W bosons, producing four W bosons and two bottom 
quarks. This analysis is sensitive to the subset of such 
events in which at least three of the W bosons decay lep- 
tonically. The heavy quark is assumed to have a mass 
of 500 GeV, corresponding to the approximate expected 
experimental limit. The normalization for this sample is 
provided at approximately next-to-next-to-leading-order 
accuracy by hathor 1.2 [42| . Both the doubly-charged 
Higgs boson and fourth-generation quark samples are 
generated with Pythia. 

The parton distribution functions for the Sherpa and 
POWHEG-BOX samples are taken from CT10 [43j], and 
from MRST2007 LO** [44] for the Pythia and her- 
wig samples. The MadGraph and ALPGEN samples use 
CTEQ6L1 [H[. The mc@nlo sample uses CTEQ6.6 [H. 

Additional pp interactions (pileup) in the same or 
nearby bunch crossings are modeled with Pythia. Simu- 
lated events are reweighted to reproduce the distribution 
of pp interactions per crossing observed in data over the 
course of the 2011 run. The mean number of interactions 
per bunch crossing for the data was ten. The luminosity 
has been measured with an uncertainty of ±3.9% (47| . 
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IV. EVENT SELECTION 

Events are required to have fired at least one single- 
electron or single-muon trigger. The electron trigger re- 
quires a minimum threshold on the momentum trans- 
verse to the beamline (px) of 20 GeV for data collected 
in the early part of 2011, and 22 GeV for data collected 
later in the year. The muon px threshold is 18 GeV for 
the full data set. The efficiency of the trigger require- 
ments for events satisfying all selection criteria ranges 
from 95% to 99% depending on the signal region, and is 
evaluated with simulated WZ events. In order to ensure 
that the efficiency is independent of the px of the lep- 
tons, the offline event selection requires that at least one 
lepton (electron or muon) has px > 25 GeV. At least one 
such lepton must also be consistent with having fired the 
relevant single-lepton trigger. A muon associated with 
the trigger must lie within \r}\ < 2.4 due to the lim- 
ited acceptance of the muon trigger, while triggered elec- 
trons must lie within < 2.47, excluding the calorime- 
ter barrel/end-cap transition region (1.37 < |?7| < 1.52). 
Additional muons in the event must lie within < 2.5 
and have px > 10 GeV. Additional electrons must sat- 
isfy the same rj requirements as triggered electrons, and 
must have px > 10 GeV. The third lepton in the event 
may be an additional electron or muon satisfying the 
same requirements as the second lepton, or a Thad with 
p T is > 15 GeV and |7? vis | < 2.5, where p T is and ?7 vis denote 
the pt and rj of the visible products of the tau decay, with 
no corrections for the momentum carried by neutrinos. 
Throughout this paper the four-momenta of tau candi- 
dates are defined only by the visible decay products. 

All parts of the detector are required to have been 
operating properly for the events under study. Events 
must have a reconstructed primary vertex candidate with 
at least three associated tracks, where each track must 
have pt > 0.4 GeV. In events with multiple primary 
vertex candidates, the primary vertex is chosen to be 
the one with the largest value of £px, where the sum is 
taken over all reconstructed tracks associated with the 
vertex. Events with pairs of leptons that are of the same 
flavor but opposite sign and have an invariant mass below 
20 GeV are excluded to avoid contributions from low- 
mass hadronic resonances. 

The lepton selection includes requirements to reduce 
the contributions from non-prompt or fake lepton candi- 
dates. These requirements exploit the transverse and lon- 
gitudinal impact parameters of their tracks with respect 
to the primary vertex, the isolation of the lepton candi- 
dates from nearby hadronic activity, and in the case of 
electron and Thad candidates, the lateral and longitudinal 
profiles of the shower in the electromagnetic calorimeter. 
There are also requirements for electrons on the quality 
of the reconstructed track and its match to the cluster 
in the calorimeter. These requirements are described in 
more detail below. 

Electron candidates are required to satisfy the "tight" 
identification criteria described in Ref. (30j . updated for 



the increased pileup in the 2011 data set. Muons must 
have tracks with hits in both the inner tracking detec- 
tor and muon spectrometer, and must satisfy criteria on 
track quality described in Ref. [31j . 

The transverse impact parameter significance is de- 
fined as \do/a(do)\, where do is the transverse impact 
parameter of the reconstructed track with respect to 
the primary vertex and cr(do) is the estimated uncer- 
tainty on do. This quantity must be less than 3.0 for 
muon candidates. Electrons must satisfy a looser cut of 
\do/<j(do)\ < 10, since interactions with material in the 
inner tracking detector often reduce the quality of the 
reconstructed track. The longitudinal impact parameter 
zq must satisfy |zosin(0)| < 1 mm for both electrons and 
muons. 

Electrons and muons are required to be isolated 
through the use of two variables sensitive to the amount 
of hadronic activity near the candidate. The first, 
PxtrackJ ^ s ^ ne scalar sum of the transverse momenta of 
all tracks with px > 1 GeV in a cone of Ai? < 0.3 around 
the lepton axis. The sum excludes the track associated 
with the lepton candidate, and also excludes tracks in- 
consistent with originating from the primary vertex. The 
second, E^° cal , is the sum of the transverse energies of 
cells in the electromagnetic and hadronic calorimeters in 
a cone of the same size. For electron candidates this sum 
excludes a rectangular region around the candidate axis 
of 0.125 x 0.172 m r] x (f> (corresponding to 5 x 7 cells in 
the main sampling layer of the electromagnetic calorime- 
ter) and is corrected for the imperfect containment of 
the electron transverse energy within the excluded re- 
gion. For muons, the sum only includes cells above a cer- 
tain threshold in order to suppress noise, and does not 
include cells with energy deposits from the muon candi- 
date. For both electrons and muons, the value of E^° cgl is 
corrected for the expected effects of pileup interactions. 
Muon candidates are required to have Px°rack/-PT < 0.13 
and -Ex°ai/PT < 0.14, while electron candidates are re- 
quired to have Px°track/.PT < 0.15 and £x° cal /px < °-14; 
see Ref. [48| for the optimization of these requirements. 

Jets in the event are reconstructed using the Fast- 
Jet [4!| implementation of the anti-fc t algorithm [Hcj . 
with distance parameter R = 0.4. The jet four-momenta 
are corrected for the non-compensating nature of the 
calorimeter, for inactive material in front of the calorime- 
ters, and for pileup |5lLl52j. Jets used in this analysis are 
required to have px > 25 GeV and lie within \rj\ < 4.9. 
Jets within the acceptance of the inner tracking detector 
must fulfill a requirement, based on tracking informa- 
tion, that they originate from the primary vertex. The 
missing transverse momentum, p™ lss , is defined as the 
negative vector sum of the transverse momenta of re- 
constructed jets, leptons, and any remaining calorimeter 
clusters unassociated with reconstructed objects. The 
magnitude of pJJ? lss is denoted E™ lss . 

Tau leptons decaying to an electron (muon) and neu- 
trinos are selected with the nominal identification cri- 
teria described above, and are classified as electrons 
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(muons). Hadronically decaying tau candidates are con- 
structed from jet candidates and are then selected using 
a boosted decision tree (BDT), which is trained to distin- 
guish hadronically decaying tau leptons from quark- and 
gluon-initiated jets 53]. The BDT is trained separately 
for tau candidates with one and three charged decay 
products, referred to as "one-prong" and "three-prong" 
taus, respectively. In this analysis, only one-prong Thad 
candidates satisfying the "tight" working point criteria 
are considered. This working point is roughly 35% effi- 
cient for one-prong Thad candidates originating from W- 
boson or Z-boson decays, and has a jet rejection factor 
of roughly 300. Additional requirements to remove Thad 
candidates initiated by prompt electrons or muons are 
also imposed. 

Since lepton and jet candidates can be reconstructed as 
multiple objects, the following logic is applied to remove 
overlaps. If two electrons are separated by AR < 0.1, 
the candidate with lower px is neglected. If a jet lies 
within AR = 0.2 of an electron or Thad candidate, the 
jet is neglected, while if the separation of the jet from an 
electron candidate satisfies 0.2 < AR < 0.4, the electron 
is neglected. In addition, electrons within AR = 0.1 of 
a muon are also neglected, as are Thad candidates within 
AR = 0.2 of electron or muon candidates. Finally, muon 
candidates with a jet within AR — 0.4 are neglected. 



V. SIGNAL REGIONS 

Events satisfying all selection criteria are classified into 
four categories. Events in which at least three of the lep- 
ton candidates are electrons or muons are selected first, 
followed by events with two electrons or muons, or one 
of each, and at least one Thad candidate. These two cat- 
egories are referred to as > 3e//i and 2e//^+ > lThad 
respectively. Next, events in each of those two cate- 
gories are sub-divided by the presence or absence of a 
reconstructed Z-boson candidate, which is defined as an 
opposite-sign same-flavor pair of lepton candidates with 
a total invariant mass within ±20 GeV of the Z-boson 
mass [HI . An additional electron may also be included in 
the combination with the same-flavor opposite-sign pair 
to satisfy the invariant mass requirement, to handle cases 
where an energetic photon from final-state radiation con- 
verts in the detector and is reconstructed as a prompt 
electron. Events with a reconstructed Z-boson candi- 
date are referred to as on-Z, and those without such a 
candidate are referred to as oS-Z. The resulting four cat- 
egories are mutually exclusive, and are chosen to isolate 
the contributions from backgrounds such as jets faking 
"Hiad candidates, and events with Z bosons produced in 
association with a jet that fakes a prompt lepton. In or- 
der to remain independent of the Z+jets control region 
described in Section IVI1 the on-Z regions have a mini- 
mum E™ 1SS requirement of 20 GeV. 

Several kinematic variables are used to characterize 
the events that satisfy all selection criteria. The vari- 



TABLE I. Kinematic signal regions defined in the analy- 
sis. The on-Z regions have an additional requirement of 
E^ iBB > 20 GeV. 

Variable Lower Bounds [GeV] Additional Requirement 
H te P ton 8 Q) 10Q) 15Q) 200i 300 

E^ isB 0, 50, 75 H^ ts < 100 GeV 

E^ iss 0, 50, 75 H% tB > 100 GeV 

m cft 0, 150, 300, 500 

m efi 0, 150, 300, 500 E^ iBB > 75 GeV 



able Hr£ p ons is defined as the scalar sum of transverse 
momenta, or p T ls for Th a d candidates, of the three lead- 
ing leptons. The variable H^ ts is defined as the sum of 
transverse momenta of all selected jets in the event. The 
"effective mass", m e ff, is the scalar sum of E™ iss , H^ ts , 
and the transverse momenta of all identified leptons in 
the event. 

Subsets of selected events are defined based on kine- 
matic properties. The 7j^P tons distribution is considered 
for all events in each category. The E!f lss distribution 
is considered separately for events with H^ ts below and 
above 100 GeV, which serves to separate events produced 
through weak and strong interactions. The m e ff distribu- 
tion is considered for events with and without a require- 
ment of E™ lss > 75 GeV. Increasing lower bounds on the 
value of each kinematic variable define signal regions; the 
lower bounds are shown in Table fl] 



VI. BACKGROUND ESTIMATION 

Standard Model processes that produce events with 
three lepton candidates fall into three classes. The first 
consists of events in which prompt leptons are produced 
in the hard interaction, including the WZ, ZZ, and 
ti + W/Z processes. A second class of events includes 
Drell-Yan production in association with an energetic 7, 
which then converts in the detector to produce a single 
reconstructed electron. A third class of events arises from 
non-prompt, non-isolated, or fake lepton candidates sat- 
isfying the identification criteria described in Section IPV1 

The first class of backgrounds is dominated by W Z —> 
tvl'l 1 and ZZ UH' events. Smaller contributions 
come from tt+W bblvl'vl"v and tt+Z -» hblvl'vl"l" 
events. Contributions from triboson events, such as 
WWW — > lvl'vl"v production, are negligible. All such 
processes are modeled with the dedicated MC samples 
described in Section Mil Reconstructed leptons in the 
simulated samples are required to be consistent with the 
decay of a vector boson or tau lepton from the hard in- 
teraction. The second class of backgrounds, from Drell- 
Yan production in association with a hard photon, is also 
modeled with MC simulation. 

The class of events that includes non-prompt or fake 
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leptons, referred to here as the reducible background, is 
estimated using in-situ techniques which rely minimally 
on simulation. Such backgrounds for muons arise from 
scmi-leptonic b- or c-hadron decays, from in-flight decays 
of pions or kaons, and from energetic jets that reach the 
muon spectrometer. Electron candidates can also arise 
from misidentified hadrons or jets. Hadronically decay- 
ing taus have large backgrounds from narrow, low-track- 
multiplicity jets that mimic Thad signatures. 

Relaxed criteria are defined for each lepton flavor. 
These criteria, in combination with a requirement that 
candidates fail the nominal identification criteria, pro- 
duce samples of lepton candidates that are rich in back- 
ground with minimal contributions from misidentified 
prompt leptons. For electrons and muons, the isolation 
criteria are relaxed to accept non-isolated leptons. Elec- 
trons are also allowed to fail the "tight" electron iden- 
tification criteria, provided they satisfy the "medium" 
criteria [30j . The relaxed Thad identification loosens the 
requirement on the BDT score. 

These samples of events are used to measure the ratio 
of the number of leptons satisfying the nominal identifica- 
tion criteria to the number that fail the nominal criteria 
but satisfy the relaxed criteria. This ratio can then be ap- 
plied as a scale factor - referred to here as a "fake factor" 
- to multilepton events satisfying the relaxed criteria to 
estimate the background in signal regions. For electron 
and muon candidates, the sample used to measure the 
fake factor consists of events that pass the high-px single- 
lepton triggers described in Section ITVl Events with more 
than one selected lepton are removed from the sample to 
avoid overlap with the signal region, and to reduce the 
contamination from Drell-Yan processes. Muons must 
also fail the nominal requirement on \do/a(do)\ to fur- 
ther remove prompt contributions. Finally, events where 
the transverse mass (tot) of the electron combined with 
the E!f lss is larger than 25 GeV are also rejected to avoid 
contamination from W+jets, where mx is defined as: 

rn T = \J (Ei + ^x iss ) 2 " IPt + PT is 1 2 - (!) 

For events with muons the transverse mass requirement 
is relaxed to 40 GeV since the inversion of the \do/a(do)\ 
requirement is sufficient to remove most of the contribu- 
tions from MK-boson decays. 

For Thad candidates, a sample of 7+jet events is used 
to measure the fake factors. The production of prompt 
photons in pp collisions is dominated by the Compton 
process qg — > (77, yielding a sample of 7+jet events that 
is rich in quark-initiated jets. In the events considered 
here, an energetic photon is used to tag the event, and 
the Thad candidate is the away-side jet. The photon is 
required to have pr > 40 GeV, and satisfy the "tight" 
identification criteria [Hj]. The photon candidate is also 
required to have E™° cgl < 5 GeV. These criteria have 
been shown to yield a mostly pure sample of photon can- 
didates, with the remainder largely consisting of events 
in which a jet fragments into a leading ir° that then de- 
cays to two photons. The resulting sample suffers from 



minimal contamination from true Thad candidates, with 
the largest contribution from W{— > Th a d^)+jet events, 
where the jet is identified as a photon, contributing less 
than 1% to the total sample. 

The fake factors for all flavors are parameterized as 
functions of the j»t and \rj\ of the candidates, to account 
for changes in the composition of the nominal and relaxed 
samples in different kinematic ranges. For electrons and 
muons with pt > 100 GeV, the fake factor is computed 
from a linear extrapolation of the fake factors between 
35 GeV and 100 GeV. An additional parameterization is 
added to account for the heavy- flavor content of the event 
based on the output of the MV1 6-tagging algorithm (5|| ; 
the largest MV1 score associated with any jet in the event 
is used to parameterize the fake factors. The correlation 
of this variable with the use of the inverted \do/a(do)\ 
requirement when estimating the muon fake factors leads 
to a bias in events with large MV1 scores, causing the 
muon fake factors to be underestimated by a factor of 
two. This bias is corrected using MC simulated samples. 

Contributions from prompt leptons can bias the re- 
ducible background estimates in two ways. The first 
arises when prompt leptons populate either the tight or 
relaxed regions when deriving the fake factors. The sec- 
ond arises when prompt leptons populate the relaxed re- 
gion when applying the fake factors. In all cases, the 
effects of prompt leptons on the reducible background 
estimates are evaluated and corrected using MC simula- 
tion. 

The background estimates are tested in several control 
regions. A control region rich in events with a Z boson 
produced in association with a jet is defined to test the 
reducible background estimates. Events in this region 
have three identified lepton candidates, with the require- 
ment that a pair of opposite-sign, same-flavor leptons has 
an invariant mass within ±20 GeV of the Z-boson mass. 
The additional requirement that the E!f lss does not ex- 
ceed 20 GeV avoids overlap with the signal regions. This 
is referred to as the Z+jets region. 

A second control region, also consisting of events with 
three lepton candidates, is defined using the low-mass 
Drell-Yan events rejected by the requirement that no 
opposite-sign, same-flavor lepton pair have m(£ + £~) < 
20 GeV. This region is referred to as the low-mass Drell- 
Yan region. 

A third region is defined in order to probe the esti- 
mates of backgrounds from non-prompt and non-isolated 
sources in events rich in heavy-flavor decays. Events 
are required to have exactly two same-sign leptons and 
^miss > 4Q Q e y Events are further required to have a 
&-jet candidate selected by the MV1 tagging algorithm, 
using a working point that is 60% efficient and that has 
a light-jet mis-tag rate of less than 1% for jets with 
Pt < 100 GeV. This sample is estimated to be primar- 
ily composed of lepton+jets it events. The same-sign re- 
quirement suppresses events where both W bosons decay 
leptonically, and enhances the contributions from events 
where one lepton candidate originates from semileptonic 



6 



TABLE II. The predicted and observed number of events 
in the Z+jets, low-mass Drell-Yan, and ti control regions. 
The Z+jets and low-mass Drell-Yan regions are populated 
by trilepton events, while the ti region is composed of same- 
sign dilepton events. Statistical and systematic uncertainties 
on the expected event yields are combined as described in 
Section ETU 



Channel 


Irreducible 


Reducible 


Total 


Observed 


Z+jets 


> 3e//j, 


165 ± 26 


160 ± 50 


320 ± 60 


359 


2e/fl+ > lThad 


3.0 ± 0.6 


1480 ± 360 


1480 ± 360 


1696 


Low-mass Drell-Yan 


> 3e//j, 


55 ± 9 


34 ± 12 


89 ± 15 


101 


2e/fl+ > lThad 


0.5 ± 0.1 


91 ± 23 


92 ± 23 


96 


tt 




25 ± 4 


58 ± 23 


83 ± 23 


87 


le/fi + lr ha d 


1.9 ± 0.4 


107 ± 27 


109 ± 27 


103 



>3 e/n ATLAS 

on-Z (s= 7 TeV 

E miss < 2Q QeV | Ldt = 4 6 fb -, 



-•-2011 Data 
Reducible 
WZ/ZZ 

■ Z+jets 

tt+V 

Syst. Unc. 




60 70 80 90 100 

Non-Z Lepton m [GeV] 



FIG. 1. The tut distribution of the _E™ 1SS and the lepton not 
associated with the Z-boson candidate decay in > 3e//i events 
in the Z+jets control region. The last bin shows the integral 
of events above 90 GeV. The bottom panel shows the ratio of 
events observed in data to those expected from background 
sources for each bin. 



b decay. This region is referred to as the ti region. An 
upper limit on H^ ts of 300 GeV reduces potential con- 
tamination from new phenomena. 

Good agreement between the expected and observed 
event yields is seen in all control regions as shown in Ta- 
ble [TTJ Figure Q] shows the tot distribution of the E™ ss 
and the lepton not associated with the Z boson candi- 
date in the > 3e//z channel of the zT+jets region. Fig- 
ure [5] shows the px distribution for the third lepton can- 
didate in the -Z+jets region. The pr distribution for 
the subleading lepton (rhad candidate) in the ti region 
is shown in Fig. |3l The pt distribution for the third 
lepton in the low-mass Drell-Yan region is shown in 
Fig. [U The H^ ptons , E^ iss , and m cfi distributions are 
not shown here, but also are in good agreement in the 
control regions. The contributions from new phenomena 
in the control regions are estimated with doubly-charged 
Higgs and fourth-generation quark events. An example 
of such contamination is shown with fourth-generation 
quark events in Fig. 3(a) where the contamination is 



small. The contributions in all other control regions from 
events with pair-produced doubly-charged Higgs bosons 
or fourth-generation quarks are negligible. 



VII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties on the predicted backgrounds 
come from several sources. These uncertainties are sum- 
marized in Table IIII1 presented as ranges of relative un- 
certainties on the total expected background yields across 
all signal regions and channels. 

The backgrounds modeled with simulated samples 
have uncertainties associated with trigger efficiencies, 
lepton efficiencies, lepton momentum scales and resolu- 



TABLE III. The range of systematic uncertainties originating 
from different sources, presented as the relative uncertainty on 
the total expected background yield in all signal regions under 
study. In cases where a source of uncertainty contributes less 
than 1% of total uncertainty in any of the signal regions, the 
minimum is presented as "< 1%". 



Source of uncertainty 


Uncertainty 


Trigger efficiency 


(< 


1) 


1% 


Electron energy scale 


(< 


1) 


- 13% 


Electron energy resolution 


(< 


1) 


1% 


Electron identification 


(< 


1) 


■ 3% 


Electron non-prompt/fake backgrounds 


(< 


1) 


- 13% 


Muon momentum scale 


(< 


1) 


1% 


Muon momentum resolution 


(< 


1) 


7% 


Muon identification 


(< 


1) 


1% 


Muon non- prompt /fake backgrounds 


(< 


1) 


- 51% 


Tau energy scale 


(< 


1) 


■ 4% 


Tau identification 


(< 


1) 


■ 4% 


Tau non- prompt /fake backgrounds 


(< 


1) 


- 24% 


Jet energy scale 


(< 


1) 


6% 


Jet energy resolution 


(< 


1) 


■ 3% 


Soft E^ iBB terms 


(< 


1) 


- 14% 


Luminosity 






3.9% 


Cross-section uncertainties 


(< 


1) 


- 14% 


Statistical uncertainties 




1 


- 25% 


Total uncertainty 




11 


- 56% 
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FIG. 2. The pt distribution of the third lepton candidate in (a) > 3e/fj, events and (b) 2e/fi+ > lrhad events in the Z+jets 
control region. The last bin in the left (right) plot shows the integral of events above 100 GeV (150 GeV). The bottom panel 
shows the ratio of events observed in data to those expected from background sources for each bin. 




FIG. 3. The pt distribution of the (a) subleading lepton in 2e/ fj, events and (b) Th a d in le//x + lThad events in the ti control 
region. The expected contribution from non-Standard-Model processes is illustrated in the left figure by events with fourth- 
generation down-type quarks (6'). The contribution from b' events in the right figure is negligible. The last bin in each plot 
shows the integral of events above 100 GeV. The bottom panel shows the ratio of events observed in data to those expected 
from background sources for each bin. 



tion, and jet energy scales and resolution. The uncer- 
tainty on the i?™ lss in simulation is computed from vary- 
ing the inputs to the E™ lss calculation within their uncer- 
tainties on the energy /momentum scale and resolution, 
and is thus strongly correlated with the other uncertain- 
ties and not presented separately. Contributions to the 
£™ ss from soft activity not associated with high-px ob- 
jects are presented separately. Uncertainties on the jet 
energy scale and resolution are significant in regions re- 
quiring large values of H^ ts or m e s, and are small oth- 
erwise. 



Uncertainties on the cross sections of the different 
Standard Model processes modeled by simulation are 
also considered. The Sherpa predictions of the WZ 
and ZZ processes are cross-checked with the next-to- 
leading-order predictions from POWHEG-box in a kine- 
matic region similar to the signal regions considered in 
this search, resulting in 10% and 25% uncertainties in the 
normalization, respectively. Uncertainties from renor- 
malization and factorization scale variations, as well as 
the variation of the parton distribution functions, con- 
tribute an additional 10% and 7% respectively, taken 
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FIG. 4. The px distribution of the (a) third lepton candidate in > 3e/fi events and (b) Thad in 2e/fi+ > lrhad events in the 
low-mass Drell-Yan control region. The last bin shows the integral of events above 50 GeV (100 GeV) in the left (right) figure. 
The bottom panel shows the ratio of events observed in data to those expected from background sources for each bin. 



from Ref. [57[. The ti+W and tt + Z backgrounds carry 
a total uncertainty of 50% based on parton distribution 
function and scale variations, and on large higher-order 
corrections [Ifl H3] • The Drell-Yan samples have a total 
uncertainty of 7% [58| . 

The reducible background estimates carry large uncer- 
tainties from several sources. A 40% uncertainty is as- 
signed to the fake factors used to estimate the reducible 
electron and muon backgrounds, based on closure stud- 
ies in MC samples and cross-checks in control regions. 
For electrons and muons with jjt > 100 GeV, where the 
fake factors are extrapolated from the values at lower px, 
a 100% uncertainty is assigned. A 100% uncertainty is 
also assigned to the fake factors for muons with high b- 
tagging scores, due to the large correction taken from MC 
simulation to remove the bias between the 6-tagging algo- 
rithm and the inverted do requirement. For the Thad fake 
estimates, a 25% uncertainty on the fake factors is deter- 
mined by altering the composition of the relaxed sample. 
In signal regions where the relaxed samples are poorly 
populated, statistical uncertainties on the reducible back- 
ground estimates become significant, especially in regions 
with high _E™ 1SS or H^! ts requirements. 

In all of the signal regions under study, the dominant 
systematic uncertainties on the total background esti- 
mate arise from the uncertainty associated with the re- 
ducible background estimates or from the uncertainty on 
the cross sections used for backgrounds taken from MC 
simulation. 

Uncertainties on the efficiency for potential sources of 
new phenomena include contributions from lepton trig- 
ger and identification efficiencies, and lepton momen- 
tum scale and resolution. Larger uncertainties on the 
signal efficiency are assigned based on variations ob- 
served between several simulated samples, including pair- 



TABLE IV. The expected and observed event yields for all 
inclusive signal channels. The expected yields are presented 
with two uncertainties, the first is the statistical uncertainty, 
and the second is the systematic uncertainty. 



Flavor Chan. 


Z Chan. 


Expected 


Observed 


> 3e//x 


off-Z 


107 ± 7 ± 24 


99 


> 3e/n 


on-Z 


510 ± 10 ± 70 


588 


2e/fl+ > lThad 


off-Z 


220 ± 5 ± 50 


226 


2e/^+ > lT h ad 


on-Z 


1060 ± 10 ± 260 


914 



production of doubly-charged Higgs bosons and of fourth- 
generation down-type quarks, and are 10% for the > 3e/fi 
channels, and 20% for the 2e//i+ > lrhad channels. 



VIII. RESULTS 

Event yields for the most inclusive signal regions in 
each search channel are presented in Table HVl No sig- 
nificant deviation from the expected background is ob- 
served. The yields for all signal regions are presented in 
Tables [VlllHXTTl of Appendix EJ 

The #,^P tons distributions for the two oS-Z signal 
channels are shown in Fig. [5l and the E™ lsa distributions 
for the same channels are shown in Fig. [6l The m e ff dis- 
tributions for the two on-Z channels are shown in Fig. [71 
The m c ff distribution for the on-Z, > 3e//i channel in 



Fig. 7(a) has 4 events with m e s > 1 TeV where a total 
of 2.2 events are expected. 

The observed event yields in different signal regions 
are used to constrain contributions from new phenomena. 
The 95% confidence level (CL) upper limits on the num- 
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FIG. 5. The H^ ptonB distribution for the off-Z (a) > 3e//i and (b) 2e/fi+ > lTbad signal channels. The dashed lines represent 
the expected contributions from events with pair-produced doubly-charged Higgs bosons with masses of 300 GeV. The last bin 
in the left (right) figure shows the integral of events above 600 GeV (500 GeV). The bottom panel shows the ratio of events 
observed in data to those expected from background sources for each bin. 




FIG. 6. The E™ lss distribution for the off-Z (a) > 3e//x and (b) 2e/p+ > lrhad signal channels. The dashed lines represent the 
expected contributions from events with fourth-generation down-type quarks with masses of 500 GeV. The last bin in the left 
(right) figure shows the integral of events above 300 GeV (200 GeV). The bottom panel shows the ratio of events observed in 
data to those expected from background sources for each bin. 



ber of events from non-Standard-Model sources (Ags) are 
calculated using the CL S method 59] . All statistical and 
systematic uncertainties on estimated backgrounds are 
incorporated into the limit-setting procedure, with corre- 
lations taken into account where appropriate. Systematic 
uncertainties on the signal efficiency are also included as 
described in Section IVIII The Ng$ limits are then con- 
verted into limits on the "visible cross section" (ogif) us- 
ing the relationship agf = -/V95/ J Ldt. 



with the median expected limits with ±la and ±2cr un- 
certainties. Observed and expected limits are also pre- 
sented in Tables ElnHXVni of Appendix [Bj The most in- 
clusive signal regions for the ff^P tons anc ] moff variables 
are composed of the same events within each channel, 
leading to identical limits. 



Figures I5HT121 show the resulting observed limits, along 
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FIG. 7. The m e ff distribution for the on-Z (a) > 3e//x and (b) 2e//i+ > lThad signal channels. The dashed lines represent the 
expected contributions from events with fourth-generation down-type quarks with masses of 500 GeV. The last bin in the left 
(right) figure shows the integral of events above 1.2 TeV (1 TeV). In the > 3e//i channel, a total of 2.13 events are expected 
for m fr > 1 TeV, and 4 events are observed. The bottom panel shows the ratio of events observed in data to those expected 
from background sources for each bin. 
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FIG. 8. The observed- and median-expected 95% CL limit 
on the visible cross section (erg") in the different signal chan- 
nels, as functions of increasing lower bounds on H^ ptonB . The 
±l<j and ±2<r uncertainties on the median expected limit are 
indicated by green and yellow bands, respectively. 
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FIG. 9. The observed- and median-expected 95% CL limit on 
the visible cross section (erg") in the different signal channels, 
as functions of increasing lower bounds on _E™ 1SS , for events 
with ffij? ts < 100 GeV. The lowest bin boundary X is GeV 
for the oS-Z channels, and 20 GeV for the on-Z channels. The 
±l<j and ±2<r uncertainties on the median expected limit are 
indicated by green and yellow bands, respectively. 



IX. MODEL TESTING 



The <7gg s limits can be converted into upper limits on 
the cross section of a specific model as follows: 

• Events from the new model are examined at the 
particle (MC-generator) level and kinematic re- 
quirements on the particles are applied. These in- 
clude the and rj requirements for leptons and 
jets, and isolation requirements for the leptons. No 



special treatment for pileup is necessary. 

• The number of events passing this selection deter- 
mines the cross section for the model given the fidu- 
cial constraints, cr fid . 

• A correction factor must be applied to take into 
account detector effects. This correction factor, 
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FIG. 10. The observed- and median-expected 95% CL limit 
on the visible cross section (agf) in the different signal chan- 
nels, as functions of increasing lower bounds on £"™ lss , for 
events with H^ s > 100 GeV. The lowest bin boundary X is 
GeV for the oS-Z channels, and 20 GeV for the on-Z chan- 
nels. The ±1<7 and ±2cr uncertainties on the median expected 
limit are indicated by green and yellow bands, respectively. 
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FIG. 12. The observed- and median-expected 95% CL limit 
on the visible cross section (a gf) in the different signal chan- 
nels, as functions of increasing lower bounds on m c g , for 
events with E™ lss > 75 GeV. The ±la and ±2<j uncertain- 
ties on the median expected limit are indicated by green and 
yellow bands, respectively. 
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FIG. 11. The observed- and median-expected 95% CL limit 
on the visible cross section (ergs 5 ) in the different signal chan- 
nels, as functions of increasing lower bounds on m c ff. The 
±1<7 and ±2<r uncertainties on the median expected limit are 
indicated by green and yellow bands, respectively. 



called e , is model-dependent, and is subject to 
uncertainties from detector resolution, reconstruc- 
tion efficiency, pileup, and vertex selection. This 
correction factor represents the ratio of the num- 
ber of events satisfying the selection criteria after 
reconstruction to all those satisfying the fiducial 
acceptance criteria at the particle level. 

• A 95% CL upper-limit on the cross section in the 
new model is then given by: 



AT ~V1S 

fid _ JV 95 _ cr 95 

The value of e fid in the > 3e//z channels ranges from 
roughly 0.50 for fourth-generation quark models to over 
0.70 for doubly-charged Higgs models producing up to 
four high-px leptons. In the 2e/fi+ > lThad channels, e fid 
is roughly 0.10 for a variety of models. Finite momentum 
resolution in the detector can cause particles with true 
momenta outside the kinematic acceptance {e.g. muons 
with pt < 10 GeV) to be accepted after reconstruction. 
The fraction of such events after selection is at most 3% 
for the > 3e//x channels and 4% for the 2e/fi+ > lrh a d 
channels. 

In order to determine e for unexplored models of new 
phenomena producing at least three prompt, isolated, 
and charged leptons in the final state, per-lepton efficien- 
cies parametrized by the lepton kinematics are provided 
here. While the experimental results are based on recon- 
structed quantities, all requirements in the following are 
defined at the particle level. 

Electrons at the particle level are required to have pt > 
10 GeV, and to satisfy |ry| < 2.47 and \r)\ £ (1.37, 1.52). 
Particle-level muons are required to have pt > 10 GeV, 
and to have \rj\ < 2.5. Electrons and muons are both 
required to be prompt, and not associated with a sec- 
ondary vertex, unless they are the product of tau-lepton 
decays. Leptonically decaying tau candidates are re- 
quired to produce electrons or muons that satisfy the 
criteria above. Hadronically decaying tau candidates are 
required to have > 10 GeV and |?? VIS | < 2.5, where 
the visible products of the tau decay include all particles 
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except neutrinos. As with reconstructed tau candidates, 
the tau four-momentum at the particle level is defined 
only by the visible decay products. 

Generated electrons and muons are further required to 
be isolated. A track isolation energy at the particle level 
corresponding to PT°track> denoted PT° truc , ^ s defined as 
the scalar sum of transverse momenta of charged particles 
within a cone of AR < 0.3 around the lepton axis. Parti- 
cles used in the sum are included after hadronization and 
must have pt > 1 GeV. A fiducial isolation energy corre- 
sponding to E^° ca \, denoted £|° trac , is defined as the sum 
of all particles inside the annulus 0.1 < AR, < 0.3 around 
the lepton axis. Neutrinos and other stable, weakly- 
interacting particles are excluded from both Px°true an( i 

Electrons must 

satisfy p?° true /p 



TABLE V. The fiducial efficiency for electrons and taus in 
different ranges pt- For tau candidates, pt = Pt s - 



3t < 0.13 and 



true ' 

/p T < 0.2, while 



muons must satisfy pf°truc/PT < 0.15 and E^ t ruc /PT < 
0.2. 

Events with at least three leptons as defined above 
must have at least two electrons and/or muons, at least 
one of which has px > 25 GeV. The third lepton is al- 
lowed to be an electron or muon, in which case the event 
is classified as a > 3e/p event, or a hadronically decaying 
tau lepton, in which case it is a 2e//i+ > lT na d event. 

A simulated sample of WZ events is used to deter- 
mine the per-lepton efficiencies eg. The leptons above are 
matched to reconstructed lepton candidates that satisfy 
the selection criteria defined in Section IIV1 with eg de- 
fined as the ratio of the number of reconstructed leptons 
satisfying all selection criteria to the number of generated 
leptons satisfying the fiducial criteria. Separate values of 
eg are measured for each lepton flavor. In the case of elec- 
trons and muons, eg is determined separately for leptons 
from tau decays. 

All efficiencies are measured as functions of the lep- 
ton pt and r\. The efficiencies for electrons and taus 
are shown in Tables [V] and IVII The 77 dependence of 
the muon efficiencies is treated by separate px efficiency 
measurements for muons with |ry| < 0.1 and those with 
|ry| > 0.1, and is shown in Table IVII1 For taus, the 
efficiency tables include the efficiency for taus gener- 
ated with p T ls < 15 GeV but reconstructed with p T ls > 
15 GeV, due to resolution effects. The corresponding effi- 
ciencies for electrons and muons generated below 10 GeV 
are much smaller, and are not included here. The final 
per-lepton efficiency for electrons and taus is obtained as 
eg = e(pr) ■ e(??)/( e )j where (e) is 0.69 for prompt elec- 
trons, 0.53 for electrons from tau decays, and 0.17 for 
hadronically decaying taus. 

The resulting per-lepton efficiencies are then combined 
to yield a selection efficiency for a given event satisfying 
the fiducial acceptance criteria. For events with exactly 
three leptons, the total efficiency for the event is the prod- 
uct of the individual lepton efficiencies. For events with 
more than three leptons, the additional leptons in order 
of descending pt only contribute to the total efficiency 
when a lepton with higher px is not selected, leading to 
terms like 616264(1 — 63), where denotes the fiducial 
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.07 
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TABLE VI. The fiducial efficiency for electrons and taus in 
different r\ ranges. For tau candidates, r\ = rj vls . 



Prompt e 



Thad 



0.0-0.1 0.675±0.003 0.52±0.01 0.210±0.009 

0.1-0.5 0.757±0.001 0.595±0.005 0.195±0.004 

0.5-1.0 0.747±0.001 0.581±0.005 0.179±0.004 

1.0-1.5 0.666±0.002 0.494±0.006 0.138±0.004 

1.5-2.0 0.607±0.002 0.465±0.006 0.170±0.004 

2.0-2.5 0.591±0.002 0.475±0.007 0.163±0.005 



efficiency for the i th pT-ordered lepton. The method can 
be extended to cover the number of leptons expected by 
the model under consideration. 

Jets at the particle level are reconstructed from all sta- 
ble particles, excluding muons and neutrinos, with the 



TABLE VII. The fiducial efficiency for muons in different px 



ranges. 


Pt 
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T — > 








[GeV] 




\v\ > 0.1 


M < 0.1 


M > 0.1 
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anti-fci algorithm using a distance parameter R = 0.4. 
Overlaps between jets and leptons are removed as de- 
scribed in Section ITVl E™ lss is denned as the magnitude 
of the vector sum of the transverse momenta of all stable, 
weakly-interacting particles, including those produced in 
models of new phenomena. The kinematic variables used 
for limit setting are defined as before: H l ° ptons i s the 
scalar sum of the transverse momenta, or p™ for r^ad 
candidates, of the three leptons that define the event; 
H^ ts is the scalar sum of all jets surviving overlap re- 
moval; _E™ 1SS is as defined above, and m e g is the sum 
of E™ lss , iJ^ ts , and all transverse momenta of selected 
leptons in the event. 

Predictions of the rate and kinematic properties of 
events with multiple leptons made with the method de- 
scribed above agree well with the same quantities after 
detector simulation for a variety of models of new phe- 
nomena. Uncertainties, based on the level of agreement 
seen across a variety of models, are estimated at 10% for 
the > 3e//x channels, and 20% for the 2e/fi+ > lri ia d 
channels. These uncertainties are included in the limits 
presented in Section [Villi 

As an example of the application of the method de- 
scribed in this section, the <jgf limits can be used to con- 
strain models predicting the pair-production of doubly- 
charged Higgs bosons. The constraints from dedicated 
and optimized analyses by ATLAS [60| and CMS [|l[ 
are expected to be stronger than the constraints obtained 
here, but these numbers serve to benchmark the results 
presented in this paper. 

Assuming a branching ratio of 100% for the decay 
i7 ±± — > ^ ± A* ± , the acceptance of the fiducial selection 
is 91% and e fid is 71% for m( J ff ±± ) = 300 GeV. The 
resulting 95% CL upper limit on the cross section times 
branching ratio (a ■ BR) is 2.5 fb. The observed and 



median expected upper limits are shown in Fig. 13(a) 



along with the observed upper limit from the dedicated 
search by ATLAS 60]. These results are obtained us- 
ing the ff}^ ptons > 300 GeV signal region in the > 3e/^, 
off-Z channel. The theoretical cross section for H ±± 
coupling to left-handed fermions (H^ ) implies that 
Hf^ masses below 330 GeV are excluded at 95% CL for 



BR(ff 



±± 



=100%. 



For the case with BR( J ff ±± ->■ ^ ± r ± ) = 100%, the ac- 
ceptance for the > 3e/fi (2e/fi+ > lrhad) channel is 24% 
(49%), and e fid is 59% (13%) for mCff^) = 200 GeV. 
The corresponding upper limit on the cross section is 12 
(19) fb, with m{H^) < 237 GeV (220 GeV) excluded 
at 95% CL. In this case, the off-Z i^?P tons > 300 GeV 
signal region is used to calculate the expected limits for 
all H ±zt masses except for m(iJ ±± )=100 GeV where the 
off-Z, H^ ptons > 200 GeV signal region is used. The 
observed and median expected limits from the > 3e//i 
channel are shown in Fig. |13(b)| 



X. CONCLUSION 

A generic search for new phenomena in events with at 
least three energetic, charged, prompt, and isolated lep- 
tons has been presented, using a data sample correspond- 
ing to an integrated luminosity of 4.6 fb -1 of pp collision 
data collected by the ATLAS experiment. The search 
was conducted in separate channels based on the pres- 
ence or absence of a hadronically decaying tau lepton or 
reconstructed Z boson, and yielded no significant devia- 
tion from background yields expected from the Standard 
Model. Upper limits at 95% confidence level on event 
yields due to non-Standard-Model processes were placed 
as a function of lower bounds on several kinematic vari- 
ables. Additional information on the fiducial selection 
of events populating the signal regions under study has 
been provided. The use of this information in the inter- 
pretation of the results in the context of models of new 
phenomena has been illustrated by setting upper limits 
on the production of doubly-charged Higgs bosons decay- 
ing to same-sign lepton pairs. 
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FIG. 13. The expected and observed 95% confidence level upper limits on the cross section times branching ratio of the (a 
H ±:iz — > fi ^ and (b) H — > fj, t final states as a function of the H ±± mass for the > 3e//i channel. For H 
the median expected limit on the ffr mass is 319 GeV and the corresponding observed limit is 330 GeV; for H 



= = 1 M ±r± i 



the median expected limit is 229 GeV and the corresponding observed limit is 237 GeV. Results from the dedicated ATLAS 
search for H — > fj, /i* 60] are also shown. 
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Appendix A: Tables of expected and observed event 
yields 



The expected and observed event yields for all signal 
regions under study are shown in Tables IVIIHIXlil 
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TABLE VIII. Results for the ^ cptons signal regions. Irreducible sources include all backgrounds estimated with MC simulation. 
Results are presented in number of expected events as N ± (statistical uncertainty) ± (systematic uncertainty). 



j^-lcptons 


Irreducible 


Reducible 


Total Exp. 


Observed 


> 3e/^, oS-Z 


GeV 


54 ±4 ±7 


54 ±6 ±23 


107 ±7 ±24 


99 


100 GeV 


32 ±2 ±4 


32 ±4 ±16 


65 ±4 ±16 


62 


150 GeV 


22 ±1 ±3 


15 ±2 ±8 


37 ±3 ±8 


27 


200 GeV 


9.7 ±0.6 ±1.5 


6 ±2 ±4 


16 ±2 ±4 


15 


300 GeV 


3.6 ±0.5 ±0.5 


2.5 ±1.2 ±1.8 


6.2 ±1.3 ±1.9 


4 


2e/M+ > lr had , oS-Z 


GeV 


6.4 ±0.4 ±1.0 


214 ±5 ±50 


220 ±5 ±50 


226 


100 GeV 


4.4 ±0.3 ±0.6 


109 ±3 ±26 


113 ±3 ±26 


113 


150 GeV 


1.7 ±0.2 ±0.3 


46 ±2 ±11 


47 ±2 ±11 


42 


200 GeV 


0.8 ±0.1 ±0.1 


17 ±1 ±4 


17 ±1 ±4 


15 


300 GeV 


0.2 ±0.1 ±0.0 


2.5 ±0.4 ±0.6 


2.7 ±0.4 ±0.6 


1 


> 3e/fi, on-Z 


GeV 


389 ±5 ±50 


120 ±8 ±40 


508 ±10 ±70 


588 


100 GeV 


285 ±4 ±40 


71 ±6 ±26 


356 ±7 ±50 


422 


150 GeV 


122 ±2 ±17 


14 ±3 ±7 


136 ±4 ±18 


151 


200 GeV 


49 ±1 ±7 


5 ±2 ±4 


54 ±2 ±8 


60 


300 GeV 


12.3 ±0.7 ±1.6 


0.5 ±0.5 ±0.5 


12.7 ±0.9 ±1.7 


18 




2e/fi+ > l-rh ad , on- 


Z 




GeV 


13.2 ±0.5 ±2.2 


1050 ±10 ±260 


1060 ±10 ±260 


914 


100 GeV 


11.1 ±0.5 ±1.9 


670 ±10 ±160 


680 ±10 ±160 


587 


150 GeV 


4.5 ±0.3 ±0.8 


66 ±2 ±16 


71 ±2 ±16 


75 


200 GeV 


1.8 ±0.2 ±0.3 


19 ±1 ±5 


21 ±1 ±5 


24 


300 GeV 


0.5 ±0.1 ±0.1 


3.0 ±0.5 ±0.8 


3.5 ±0.5 ±0.8 


7 



TABLE IX. Results for the _E™ 1SS , _H^T ts <100 GeV signal regions. Irreducible sources include all backgrounds estimated with 
MC simulation. Results are presented in number of expected events as A?" ± (statistical uncertainty) ± (systematic uncertainty). 



E™ iBS > 


Irreducible 


Reducible 


Total Exp. 


Observed 


> Se/n, oS-Z 


GeV 


46 ±4 ±6 


41 ±5 ±16 


86 ±6 ±17 


89 


20 GeV 


28 ±4 ±3 


28 ±4 ±12 


56 ±6 ±12 


65 


50 GeV 


7.5 ±0.5 ±1.0 


15 ±2 ±7 


22 ±2 ±7 


25 


75 GeV 


3.0 ±0.3 ±0.4 


7 ±2 ±4 


10 ±2 ±4 


10 


2e/M± > lr had , oS-Z 


GeV 


5.3 ±0.4 ±0.9 


184 ±4 ±40 


190 ±4 ±40 


202 


20 GeV 


4.4 ±0.3 ±0.7 


93 ±3 ±20 


98 ±3 ±20 


91 


50 GeV 


1.5 ±0.2 ±0.2 


17 ±1 ±4 


19 ±1 ±4 


20 


75 GeV 


0.6 ±0.1 ±0.1 


8.0 ±0.8 ±1.8 


8.5 ±0.8 ±1.8 


10 


> 3e/fi, on-Z 


20 GeV 


340 ±5 ±50 


100 ±7 ±31 


439 ±9 ±60 


509 


50 GeV 


105 ±2 ±14 


14 ±3 ±5 


119 ±3 ±14 


144 


75 GeV 


40 ±1 ±5 


5 ±1 ±2 


46 ±2 ±6 


57 


2e/fi+ > l-r had , on-Z 


20 GeV 


11.3 ±0.5 ±1.9 


984 ±10 ±240 


1000 ±10 ±240 


862 


50 GeV 


4.6 ±0.3 ±0.7 


43 ±2 ±11 


48 ±2 ±11 


33 


75 GeV 


2.0 ±0.2 ±0.3 


4.1 ±0.6 ±1.0 


6.1 ±0.6 ±1.0 


4 
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TABLE X. Results for the E^ iBB , H^ ts > lOOGeV signal regions. Irreducible sources include all backgrounds estimated with MC 
simulation. Results are presented in number of expected events as N ± (statistical uncertainty) ± (systematic uncertainty). 



E™ iss > 


Irreducible 


Reducible 


Total Exp. 


Observed 


> 3e//i, oS-Z 


GeV 


7.7 ±0.8 ±1.2 


13 ±2 ±7 


21 ±2 ±7 


10 


20 GeV 


6.0 ±0.6 ±0.9 


12 ±2 ±6 


18 ±2 ±6 


8 


50 GeV 


3.2 ±0.3 ±0.5 


8 ±2 ±5 


11 ±2 ±5 


5 


75 GeV 


2.2 ±0.2 ±0.3 


7 ±2 ±4 


9 ±2 ±4 


5 


2e/M+ > iThad, oS-Z 


GeV 


1.1 ±0.1 ±0.2 


30 ±2 ±7 


31 ±2 ±7 


24 


20 GeV 


1.1 ±0.1 ±0.2 


23 ±1 ±6 


25 ±1 ±6 


20 


50 GeV 


0.7 ±0.1 ±0.1 


14.5 ±1.1 ±3.4 


15.2 ±1.1 ±3.4 


13 


75 GeV 


0.5 ±0.1 ±0.1 


9.3 ±0.8 ±2.2 


9.8 ±0.8 ±2.3 


8 


> 3e/fi, on-Z 


20 GeV 


49 ±1 ±7 


20 ±4 ±10 


69 ±4 ±12 


79 


50 GeV 


29 ±1 ±4 


7 ±2 ±3 


36 ±2 ±5 


43 


75 GeV 


17.4 ±0.7 ±2.1 


5 ±1 ±2 


22 ±2 ±3 


28 


2e/|U+ > lTdad, on-Z 


20 GeV 


1.9 ±0.2 ±0.4 


61 ±2 ±15 


63 ±2 ±15 


52 


50 GeV 


1.1 ±0.1 ±0.2 


7.8 ±0.8 ±1.9 


8.9 ±0.8 ±1.9 


11 


75 GeV 


0.7 ±0.1 ±0.1 


2.7 ±0.4 ±0.7 


3.4 ±0.5 ±0.7 


1 



TABLE XL Results for the m e s signal regions. Irreducible sources include all backgrounds estimated with MC simulation. 
Results are presented in number of expected events as N ± (statistical uncertainty) ± (systematic uncertainty). 

m e ff > Irreducible Reducible Total Exp. Observed 



> 3e/n, oS-Z 



GeV 54 ±4 ±7 54 ±6 ±23 107 ±7 ±24 99 

150 GeV 32 ±2 ±4 43 ±4 ±20 75 ±4 ±20 64 

300 GeV 12.0 ±0.9 ±1.6 16 ±2 ±8 28 ±3 ±8 15 

500 GeV 3.3 ±0.2 ±0.5 3.2 ±1.2 ±2.4 6.5 ±1.2 ±2.5 5 



2e/M± > 1-rhad, oS-Z 



GeV 6.4 ±0.4 ±1.0 214 ±5 ±50 220 ±5 ±50 226 

150 GeV 4.4 ±0.3 ±0.7 106 ±3 ±24 111 ±3 ±24 101 

300 GeV 1.3 ±0.2 ±0.2 31 ±2 ±7 32 ±2 ±7 25 

500 GeV 0.4 ±0.1 ±0.2 6.6 ±0.7 ±1.6 7.0 ±0.7 ±1.6 6 



> 3e/fi, on-Z 



GeV 390 ±5 ±50 120 ±8 ±40 510 ±10 ±70 588 

150 GeV 270 ±3 ±40 57 ±6 ±22 330 ±7 ±40 399 

300 GeV 73 ±1 ±10 16 ±3 ±8 89 ±4 ±13 103 

500 GeV 22.2 ±0.9 ±2.8 3 ±1 ±1 25 ±2 ±3 29 



2e/^+ > l7i md , on-Z 



GeV 13.2 ±0.5 ±2.2 1050 ±10 ±260 1060 ±10 ±260 914 

150 GeV 10.7 ±0.5 ±1.8 360 ±5 ±90 370 ±5 ±90 309 

300 GeV 2.9 ±0.3 ±0.4 47 ±2 ±12 50 ±2 ±12 42 

500 GeV 0.9 ±0.2 ±0.1 7.7 ±0.8 ±1.9 8.7 ±0.8 ±2.0 5 
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TABLE XII. Results for the m c ff, high-_E™ lss signal regions. Irreducible sources include all backgrounds estimated with MC 
simulation. Results are presented in number of expected events as N ± (statistical uncertainty) ± (systematic uncertainty). 

m e ff > Irreducible Reducible Total Exp. Observed 
> 3e//j, oS-Z 

GcV 5.1 ±0.4 ±0.7 13 ±2 ±8 18 ±2 ±8 15 

150 GeV 5.1 ±0.4 ±0.7 13 ±2 ±8 18 ±2 ±8 15 

300 GcV 3.7 ±0.3 ±0.5 10 ±2 ±6 13 ±2 ±6 9 

500 GeV 1.7 ±0.2 ±0.2 2.9 ±1.1 ±2.3 4.5 ±1.1 ±2.3 4 
2e/M± > lrhad, oS-Z 

GeV 1.0 ±0.2 ±0.1 17 ±1 ±4 18 ±1 ±4 18 

150 GcV 1.0 ±0.2 ±0.1 17 ±1 ±4 18 ±1 ±4 18 

300 GeV 0.6 ±0.1 ±0.1 11.9 ±0.9 ±2.9 12.4 ±0.9 ±2.9 11 

500 GcV 0.2 ±0.1 ±0.1 3.2 ±0.5 ±0.8 3.4 ±0.5 ±0.8 2 
> 3e/n, on-Z 

GeV 58 ±1 ±7 10 ±2 ±4 68 ±2 ±8 85 

150 GeV 58 ±1 ±7 10 ±2 ±4 68 ±2 ±8 85 

300 GeV 32 ±1 ±4 6 ±1 ±2 37 ±2 ±4 47 

500 GcV 11.8 ±0.6 ±1.4 2.2 ±1.1 ±0.7 14.0 ±1.3 ±1.6 18 

2e/^+ > lr had , on-Z 

GeV 2.7 ±0.3 ±0.4 6.8 ±0.7 ±1.6 9.5 ±0.8 ±1.7 5 
150 GcV 2.7 ±0.3 ±0.4 6.7 ±0.7 ±1.6 9.4 ±0.8 ±1.7 4 
300 GcV 1.6 ±0.2 ±0.2 3.5 ±0.5 ±0.9 5.0 ±0.5 ±0.9 2 
500 GeV 0.6 ±0.1 ±0.1 0.4 ±0.1 ±0.1 1.0 ±0.2 ±0.1 
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Appendix B: Tables of expected and observed limits 



The expected and observed 95% confidence level up- 
per limits on the expected event yields from new phe- 
nomena for all signal regions under study are shown in 
Tables IXTTIllXVlTl 



TABLE XIII. Limits in the H^ pton; 



bins shown as the upper 
limit on the visible cross section (o-gJf = -/V95/ J Ldt). 



j^-lcptons 

[GeV] 



Observed Expected +1(T _ iCT +2<T 
[fb] [fb] [fb] [fb] 



> 3e/n oS-Z 


> 


11 


11 


5 


9 




> 100 


8.7 


8.5 


2.9 

1.6 


6.9 


2.G 


> 150 


4.0 


4.6 


1.8 

1.2 


5.1 


1.9 


> 200 


4.4 


3.6 


1.7 

1.0 


4.9 


1.3 


> 300 


1.6 


1.9 


1.0 

0.4 


2.4 


O.G 




2e/jLt+ > 1-Thad 


off-Z 






> 


25 


23 


13 


29 


9 


> 100 


14 


14 


6 

3 


17 


5 


> 150 


6.1 


6.4 


3.4 

1.3 


8.3 


:'>. 1 


> 200 


3.3 


3.6 


1.9 

1.2 


5.0 


1.7 


> 300 


1.2 


1.5 


1.0 

0.5 


2.4 


0.8 






>3e/)i on- 


■Z 






> 


48 


33 


15 

g 


32 


14 


> 100 


38 


25 


11 

7 


23 




> 150 


14 


12 


4 

4 


10 


5 


> 200 


7.2 


6.5 


2.8 

1.9 


6.2 


2.6 


> 300 


4.5 


3.1 


1.4 

0.7 


3.5 


1.1 




2e//i+ > lThad 


on-Z 






> 


85 


94 


41 

22 


96 


30 


> 100 


53 


61 


26 

16 


64 


27 


> 150 


11.0 


9.9 


4.3 

2.2 


11.0 


3.7 


> 200 


5.2 


4.5 


2.0 

1.3 


5.3 


1.9 


> 300 


3.0 


1.9 


1.0 

0.6 


2.7 


0.8 



TABLE XIV. Limits in the E^ iss bins with H^ tB > 100 GeV 
requirement shown as the upper limit on the visible cross 
section (a^f = N 95 / J Ldt). 



rpmiss 


Observed Expected 


+ — 1(7 


+ 2(T 

-2<T 


[GeV] 


[fb] 


[fb] 


[fb] 


[fb] 


> 3e/n oS-Z 


> 


2.6 


3.1 


1.5 

0.7 


3.4 

1.4 


> 50 


2.1 


2.4 


1.0 

0.8 


2.3 

1.2 


> 75 


2.1 


2.3 


1.1 

0.4 


1.9 

0.9 


2e/n+ > lT had oS-Z 


> 


4.2 


4.8 


2.5 

1.5 


6.1 

2.1 


> 50 


3.1 


3.3 


1.8 

1.2 


4.4 

1.6 


> 75 


2.6 


2.1 


0.8 

0.6 


1.9 

1.0 






> 3e/n on- 


z 




> 20 


11.0 


8.7 


2.5 

3.2 


7.0 

4.1 


> 50 


6.4 


4.9 


2.3 

1.2 


5.4 

1.8 


> 75 


5.1 


3.8 


1.6 

1.0 


3.8 

1.4 




2e/n+ > lThad 


oa-Z 




> 20 


5.9 


7.3 


2.9 

1.4 


6.8 

3.5 


> 50 


3.4 


2.8 


1.2 

0.9 


3.8 

1.2 


> 75 


1.2 


1.5 


0.4 

0.4 


1.0 

0.6 



TABLE XV. Limits in the E^ 1SB bins with H^ ta < 100 GeV 
requirement shown as the upper limit on the visible cross 
section (a^f = N g5 / j Ldt). 



TTimiss 


Observed Expected 




+2cr 

-2ct 


[GeV] 


[fb] 


[fb] 


[fb] 


[fb] 


> 3e//x oS-Z 


> 


11 


10 


4 

2 


8 

4 


> 50 


5.3 


4.7 


1.9 

1.0 


4.5 

1.6 


> 75 


3.1 


3.0 


1.0 

0.6 


1.8 

1.0 




2e/n+ > lr had 


oft-Z 




> 


23 


21 


9 

6 


23 

9 


> 50 


4.3 


4.0 


2.3 

1.2 


5.0 

1.7 


> 75 


3.1 


2.6 


1.1 

0.7 


3.1 

1.1 






> 3e/n on- 


z 




> 20 


41 


30 


10 

9 


20 

14 


> 50 


16 


10 


4 

3 


11 

5 


> 75 


8.0 


5.4 


2.6 

1.3 


6.2 

1.9 




2e/n+ > lThad 


oa-Z 




> 20 


80 


88 


39 

23 


94 

43 


> 50 


4.4 


5.5 


3.2 

1.4 


7.6 

2.1 


> 75 


1.8 


2.2 


0.4 

0.7 


1.0 

1.0 
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TABLE XVI. Limits in the m e g bins shown as the upper limit 
on the visible cross section (crJJf = Ng$/ J Ldt). 





Observed Expected 


-la 


+ 2CT -2a 


[GeV] 


[fb] 


[fb] 


[fb] 


[fb] 


> 3e/n oS-Z 


> 


11 


11 


5 

2 


9 

4 


> 150 


8.1 


8.8 


3.0 

2.2 


7.2 

3.9 


> 300 


3.1 


3.7 


1.7 

0.7 


3.8 

1.6 


> 500 


2.1 


2.1 


1.1 

0.6 


2.3 

0.9 


2e/n+ > lT had oS-Z 


> 


25 


23 


13 

5 


29 

9 


> 150 


12 


13 


6 

4 


14 

5 


> 300 


3.9 


4.9 


2.5 

1.5 


6.4 

2.3 


> 500 


2.2 


2.4 


1.3 

0.5 


3.4 

1.2 






> 3e/n on- 


z 




> 


48 


33 


15 

8 


32 

14 


> 150 


37 


25 


9 

7 


21 

11 


> 300 


11 


9 


4 

2 


9 

3 


> 500 


4.8 


3.9 


1.7 

1.0 


4.3 

1.1 




2e/n+ > lThad 


oa-Z 




> 


85 


94 


41 

22 


96 

30 


> 150 


28 


35 


13 

11 


34 

15 


> 300 


5.9 


6.8 


2.8 

1.8 


8.1 

2.4 


> 500 


1.9 


2.5 


1.4 

1.0 


3.5 

1.2 
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requirement shown as the upper limit on the visible cross 
section (<7^ s = JV 95 / / Ldt). 



+ 1(7 +2t7 



[GeV] [fb] [fb] [fb] [fb] 

> 3e/fi oS-Z 



> 


3.8 


3.9 


1.5 

0.7 


3.4 


1.3 


> 150 


3.8 


3.9 


1.5 

1.0 


3.6 


1.3 


> 300 


2.8 


3.0 


1.2 

0.7 


3.2 


1.1 


> 500 


2.1 


2.0 


0.8 

0.4 


2.2 


0.8 




2e//i+ > lr had 


oft-Z 






> 


3.9 


3.8 


2.1 

1.2 


5.2 


1.8 


> 150 


4.0 


3.9 


2.0 

1.3 


4.8 


1.7 


> 300 


2.9 


3.1 


1.6 

1.0 


4.1 


1.3 


> 500 


1.5 


1.6 


0.7 

0.5 


2.1 


0.7 






> 3e/n on- 


z 






> 


10.0 


6.9 


3.0 

1.3 


7.4 


2.5 


> 150 


10.0 


7.1 


2.8 

2.2 


7.0 


2.5 


> 300 


6.8 


4.9 


2.1 

1.0 


5.1 


2.1 


> 500 


3.9 


3.0 


1.2 

0.7 


3.4 


1.3 




2e/n+ > lr ha d 


on-Z 






> 


1.6 


2.4 


1.4 

0.9 


3.8 


1.3 


> 150 


1.4 


2.5 


1.5 

0.8 


3.8 


1.5 


> 300 


1.5 


2.0 


1.1 

0.8 


2.7 


1.1 


> 500 


0.9 


1.1 


0.8 

0.4 


2.0 


0.4 
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